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The concept of the half density matrix is proposed. It unifies the quantum states which are described by 
density matrices and physical processes which are described by completely positive maps. With the help of 
the half-density-matrix representation of Hermitian linear map, we show that every positive map which is not 
completely positive is a difference of two completely positive maps. A necessary and sufficient condition for a 
positive map which is not completely positive is also presented, which is illustrated by some examples. 
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Entanglement has become one of the central concept in 
quantum mechanics, specially in quantum information. A 
quantum state of a bipartite system is entangled if it cannot be 
prepared locally or it cannot be expressed as a convex combi- 
nation of direct product states of two subsystems. This kind 
of state is also called inseparable. Though easily defined, it is 
very hard to recognize the inseparability of a mixed state of a 
bipartite system. 

An operational-friendly criterion of separability was pro- 
posed by Peres [|l]]. This criterion is based on the observa- 
tion that the partial transposition of a separable density matrix 
remains positive semidefinite. That the partial transposition 
of a density matrix is not positive semidefinite infers the in- 
separability of the density matrix. This provides a necessary 
condition for the separability. There exist entangled states 
with positive partial transposition, which exhibit bound en- 
tanglement [^]. Examples of such kind were first provided in 
Ref. ^ and then constructed in Ref. systematically with 
the help of unextendible product basis. 

Later on, by noticing that the transposition is a positive map 
(to be described later in details), a necessary and sufficient 
condition of the separability was proposed in Ref. [^J: A bi- 
partite state is separable iff it is still positive semidefinite un- 
der all positive maps acting on a subsystem. In other words, a 
density matrix of a bipartite system is inseparable iff there ex- 
ists a positive map acting on a subsystem such that the image 
of the density matrix is not positive semidefinite. Hence the 
inseparability can be recognized by positive maps which are 
not completely positive. 

Completely positive maps, which are able to describe the 
most general physical process [^, are better understood than 
positive maps which are not completely positive. Positive 
maps from Hilbert space 7^2 (two-dimensional) to 7^2 or Tij, 
are all decomposable [^, which are characterized by transpo- 
sition and completely positive maps only. As a result in the 
cases of 7^2 x 7^2 and 7i2 x ^3 the transposition criterion is 
also a sufficient condition for separability |^ . Therefore fur- 
ther understandings of positive maps which are not completely 
positive will facilitate the recognition and classification of the 
inseparable mixed states. 

As a direct calculation will show, under an orthonormal and 
complete basis { |n) }^Zq the transposition of an L x L matrix 



p can be expressed as 

L-l 

= Trp - ^ Gran P (T^n^ (1) 
m.n— 

where cTmn = — \n){m\)/^/2. We see immediately 

that the transposition is a difference of two completely pos- 
itive maps. And this statement will be proved to hold true 
for all positive maps which are not completely positive, which 
will be also characterized by a necessary and sufficient condi- 
tion in this Letter 

For this purpose we shall first develop an extremely useful 
tool — half density matrix that unifies the description of the 
quantum states and physical processes. And then we derive 
a half-density-matrix representation of an arbitrary Hermitian 
linear map from which our main results are obtained. Along 
with the introduction of the concept of half density matrix, its 
relations to the ensembles and the purifications of mixed states 
are clarified and its applications in the field of quantum infor- 
mation such as quantum teleportation are also presented. 

Normally, quantum states, pure or mixed, are described 
by density matrices, positive semidefinite operators (whose 
eigenvalues are all nonnegative) on the Hilbert space of the 
system. Because of its property of positive semidefinite the 
density matrix p can always be written as p = TT^ where 
matrix T is called here as the half density matrix (HDM) for a 
quantum state. 

Obviously, the half density matrix for a given density ma- 
trix is not unique. For example TU and T are correspond- 
ing to the same mixed state p = TT^ whenever U is uni- 
tary. Generally, the half density matrix T for a mixed state 
p of an s-level system is an s x L rectangular matrix with 
L > r = Rank(p), i.e., a linear map from an L-dimensional 
Hilbert space Hl to an s-dimensional Hilbert space Hs. The 
rank r of the density matrix equals to the rank of the half den- 
sity matrix T and r = 1 for pure state. 

Under an orthonormal and complete bases {|?ti)}^^o ^'^'^ 
{|n)}^^Q of Hilbert spaces Hs and 7Yl, a typical half den- 
sity matrix of dimension s x L can be constructed as = 
yt(As, Osx(L-s)) 1 where is a diagonal s x s matrix 
formed by all the square roots of the eigenvalues of p (the 
singular numbers of Tg) and V is an s x s unitary matrix diag- 
onalizing the density matrix p. Obviously we have p = TeTl 
As a direct result of the singular number decomposition of an 
arbitrary matrix we have the following 
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Lemma: Given a density matrix p of an s-level system, an 
s X L matrix T is a half density matrix for p, i.e., p ~ TT\ 
if and only if there exists an L x L unitary matrix U such that 

T = TeU. 

When written explicitly in the established basis, the rela- 
tion p = TeT^ results in exactly an ensemble formed by all 
the eigenvectors X^t|y7j^ of the mixed state, which is referred 
to as eigen-ensemble here. In this way every half density ma- 
trix T of a mixed state p corresponds to an ensemble of the 
mixed state. The above Lemma tells us that every ensemble 
of a given mixed state is related to the eigen-ensemble by a 
unitary matrix which has been proved by other means JlO|]. 
Therefore the half density matrix of a density matrix is phys- 
ically equivalent to an ensemble of the corresponding mixed 
state. 

Every mixed state p admits a purification [[TTI], a pure state 
10) of a bipartite system including this system as a subsys- 
tem such that p — Tr2|(/)) ((/>|. Under the established basis, a 
general pure state in is 

\4>) = Y.Y. C-^n\m)^\n)2 T\<^l)- (2) 

m— n— 

Here pure state | $i ) ^ ~q | i | f^) 2 lives in Hilbert space 
TCl^'Hl and T is a linear map from Hl to Hs acting on the 
first I/-dimensional Hilbert space Hl- Under the given bases 
linear map T is represented by an s x L rectangular matrix 
with matrix elements given by (m|T|n) = Cmn- When the 
pure state is normalized we have Tr(T^r) = 1. Alterna- 
tively, we also have |0) = r^|$s) with state \^s) defined in 
Hs^'Hs similar to state |"I>l). The linear mapT^ -.Hs^Hl 
acts on the second Hs and it is represented by the transposi- 
tion of T under the established basis. 

Tracing out the second system we obtain a reduced density 
matrix of the first subsystem ps = TT^ and similarly pl = 
rpTrp* Jqj. jjjg gg(.Qjj(j subsystem. That is to say T is the HDM 

for the reduced half density matrix ps = Tr/, |0) ((/>| of the first 
subsystem and its transposition T-^ for pl — Tr^lt/)) {(f)\. Thus 
a one-to-one correspondence between a normalized pure state 
10) of a bipartite system, a purification, and a linear map T 
satisfying Tr(T^r) = 1, a half density matrix, is established. 
Therefore a half density matrix T is also equivalent to a pu- 
rification of the mixed state. The linear map T is also referred 
to as the half density matrix of a bipartite pure state, which is 
unique by definition. If s = L the polar decomposition of T 
will result in the useful Schmit-decomposition. 

The pure bipartite state is separable iff the rank of its half 
density matrix is one. For a pure product state \v)s\w) l the 
half density matrix is |u)(u)*| where |w*) is the index state 
of state \w) defined by |w*) = {w\^l) [p^. For later use we 
dtfms a. mirror operator M L = \ in the Hilbert space 

Hl ® Hl, which has the property {w*\Ml\w*) — |w)(w|. 
The partial transposition of the mirror operator X = 7\/J^ is 
in fact the exchanging (or swapping) operator introduced by 
Werner [ p^ (denoted as V there). 

As an application, we consider a state \(f>) ul'fp) 3 of a tri- 
partite system with all three subsystems 1,2 and 3 being s- 
level systems. Let denote the HDM of the bipartite state 



10)12 and \k; 1)23 = Tki\^s)23 denote an orthonormal com- 
plete basis for systems 2 and 3 with HDMs Tki satisfying 

TiTkiTli, = 5kk'5w for orthogonality and TfeiOT^^; = 
TrO for completeness. We then have expansion 

|</')l2|^)3= E ^t'^M IV')l|fc;0 23. (3) 

fc,;=o 

This describes exactly a quantum teleportation of an unknown 
quantum state |-0) from system 3 to system 1 when both 
and Tki are unitary or state |0)i2 and basis |fc; 1)23 are maxi- 
mally entangled states |]l3|]. 

The mixed state psL of an (s x L) bipartite system can also 
be equivalently and conveniently characterized by HDMs of 
pure bipartite states. Let {|0i),Pi}fLi be an ensembles of psL 
we have 

-R R 
PsL =Y.P^\(i^^){<P^\ = ^ AA^i^, (4) 
i=l i=l 

where we have denoted Ai as the half density matrix of the 
pure state i.e., ^/pl|0i) = Ai|$L)- Obviously 

HDMs defined by Ai = Uij Aj characterize the same den- 
sity matrix whenever J7 is a unitary matrix. And from the 
Lemma we know that given a density matrix this is the only 
freedom that the half density matrices can have. 

The density matrix expressed in the form as in Eq.(^ can be 
easily manipulated by local operations. For example the den- 
sity matrix under operation Us ® is transformed to density 
matrix specified by half density matrices UsAiUl^. The tilde 
operation p ^ p introduced in Ref. [ |l4| ] to obtain explicitly 
the entanglement of formation of two-qubit is simply an anti- 
linear transformation Ai — TrA|^ — Al[. 

In the discussions above we have defined the half density 
matrices for the states of a single system, for pure bipartite 
states, and for mixed bipartite states. The physical processes 
can also be characterized by half density matrices. A general 
physical process which can include unitary evolutions, trac- 
ing out one system, and general measurements is described 
by trace-preserving completely positive maps ^TJ, which 
is a special kind of Hermitian linear map. 

A Hermitian linear map sends linearly Hermitian opera- 
tors to Hermitian operators that may live in different Hilbert 
spaces. Let C denote a general Hermitian linear map from 
Hilbert space Hl Hs- Because the map £ is linear the 
map £ is also a Hermitian linear map from Hl ^ Hl 
to Hs <8i Hl, where Tl denotes the identity map on Hl- Re- 
calling that the mirror operator AIl = |$l) ($l | is defined on 
Hl <8) Hl, its image 

HsL^ C®Il{Ml) (5) 

is therefore a Hermitian operator in Hs ® Hl- Let |-0j^) = 
Ai\^L) (i < i+) denote the eigenvectors corresponding to 
the positive eigenvalues of i/sL and |^/;~) = Bi\^L) {i < 
to negative eigenvalues of Hsl, where i± is the number of the 
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positive/negative eigenvalues of Hsl- We then have 

HsL - ^^^^lA\ - B^MLBl , (6) 

i=l i=l 

in which the norms of the eigenvectors \'4jf) have been taken 
to be the absolute values of corresponding eigenvalues. Be- 
cause the eigenvectors corresponding to different eigenvalues 
are orthonormal we have Tr(Aii3j) = for all i and j. In this 
sense two families of half density matrices {Ai] and {Bi\ are 
orthogonal to each other. 

For a pure state P^, — \w){w\ in the Hilbert space Hl we 
have — {w*\M]^\w*) where \w*) is the index state of \w). 
As a result we have C{Pw) — {w*\Hsl\w*) . Taking into 
account of the hnearity of the Hermitian map C, we finally 
obtain 

C{H) - " Y B^HBl (7) 

i=l i=l 

where is an arbitrary Hermitian matrix in Hl- This is called 
the half-density-matrix representation of a Hermitian linear 
map. As one result we have 

® C{^sl)\^s) = Tt{Hsl^Jl) (8) 

for an arbitrary Hermitian matrix in TCs ® Hl- As an- 
other consequence, a one-to-one correspondence between the 
Hermitian maps C : Hl '—^ Ti-s and Hermitian matrix Hsl 
(an observable) in Hs can be established 

C{H)=Tvl{HslH^) (9) 

in addition to Eq.(|]). 

The HDM representation of Hermitian linear map is not 
unique. Suppose two integers M > i+ and N > and 
let SU (M, N) denote the pseudo-unitary group formed by 
(M + A^) X {M + N) matrices satisfying SrjS^ — r\ where 
r] = ImQ) (-In) and Im(n) is the M x M (N x N) identity 
matrix. If we define a family of HDMs {T,}f£+^ as = A, 
(1 < i < «+), = Bi {M + 1 < i < M + i_) and = 
otherwise and take an arbitrary element S of SU{M, N), a 
new family of HDMs {f,}fi+^ defined by f, = J2. S^jTj 
represents the same Hermitian linear map 

M N 

C{H) = Y ~ Y ^J^^l (10) 

i=l j=l 

A positive map is a special Hermitian linear map which 
maps any positive semidefinite operator to a positive semidef- 
inite operator A Hermitian linear map S : Hl '-^ Hs is 
positive if and only ifTr(5(gL)Ps) = Tr ( i?J^ (g) Ql) > 
for all pure product state P, (g) Ql where — S®Il{Ml)- 
In the following S is always a positive map. 

A completely positive (CP) map is a positive map which 
keeps its positivity when the system it acts on is embedded as 
a subsystem in an arbitrary larger system. That is, for a CP 



map S : Hl ^ Hs and an arbitrary positive integer k the 
induced map S ®Ik from Ti/, (g) Hk to Hs ® Hk is positive. 

However it is enough to check whether the image Hsl — 
S (g Il{Ml) of the mirror operator Ml is positive semidefi- 
nite or not. If it is positive semidefinite, then the negative part 
in the HDM representation Eq.(0) disappears, which yields 
exactly the operator-sum representation of a CP map nTIn 

*+ 

S{p) = YA.pA\. (11) 

1=1 

If the trace is preserved, we have further A\Ai — 1. 
Therefore the operator-sum representation of a CP map can 
also be referred to as a half-density-matrix representation. Es- 
pecially, if HsL equals to the identity matrix Is ® II, the cor- 
responding CP map is simply the trace operation St{p) = 
/sTrp. 

A positive map which is not completely positive (non-CP) 
is nonetheless a Hermitian map so that it has a HDM represen- 
tation as Eq.(^, from which we obtain S = Sa — Sb where 
two CP maps Sa.b are represented by HDMs {Ai\ and {Bi] 
respectively. Two CP maps Sa,b are said to be orthogonal if 
their HDMs are orthogonal to each other, i.e., Tr(v4ii3j) = 
for all We see that Hsl can not be positive semidefinite. 

Conversely, if the Hermitian matrix Hsl has at least one 
negative eigenvalue then it determines a non-CP positive map. 
Let I -tp) denote an eigenvector corresponding to one of the neg- 
ative eigenvalues of Hsl and P^ = \'ip) From identity (g|) 
we see immediately thatX^ (g S{P^) is not positive semidef- 
inite, i.e., map S is not completely positive. We note that the 
eigenspace corresponding to the negative eigenvalues of Hsl 
contains no product state because of positivity. To summarize, 
we have the following 

Theorem: Every positive map which is not completely pos- 
itive is a difference of two orthogonal completely positive 
maps; A Hermitian linear map S : Hl ^ Hs is positive 
but not completely positive if and only if for all pure product 
state Ps (g Ql in Hl <» Hs we have Tx{HslPs ® Ql) > 
while HsL = 5 (g Xl{Ml) is nof positive semidefinite. 

This theorem provides an obvious way to construct a non- 
CP positive map form Hl to Hs- First, we choose a proper 
Hermitian matrix Hsl in Hs (g Hl satisfying the conditions 
specified in the above theorem. Then a non-CP positive map 
S : Hl Hs is determined by S{pl) — TrL{HsLpJ,)- 

As the first example we consider the the exchanging opera- 
tor defined in Hl g) Hl hy X = mJ^^ or explicitly 

L-l 

X= Y \m,n}{n,m\. (12) 

m,n=0 

The exchanging operator X has two eigenvalues ±1 and 
o'mnl^L) {iTi > n) are the eigenvectors corresponding to 
eigenvalue —1. Therefore X is not positive semidefinite 
and for any pure product states \pp) — \v)\w) we have 
{pp\X\pp) = \{v\w)\'^ > as specified by the above theorem. 
In fact the resulting non-CP positive map on TY^ is exactly the 
transposition = Tr2(X/9^). By writing X in its diagonal 
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form we obtain = St{p) — Sa{p), where the CP map Sa 
is represented by HDMs {(Tmn} and St is the trace operation. 

As the second example we consider a Hermitian matrix in 
TLl ^Ti-L defined by Hn — II ® II — M^. It is not positive 
semidefinite because {^l\Hii\^l) < and for every product 
states \pp) we have {pp\Hii\pp) = 1 — > 0. Ac- 

cordingly, a non-CP positive map is defined on Hl as A(p) = 
Tip — p, which provides the reduction criterion Jl^ , 16 1: Ev- 
ery inseparable state in Ti^ ® Hl which loses its positivity 
under map X^, "X" A is distillable and in the distillation pro- 
cedure provided in Ref. [|l5| the HDM of pure bipartite state 
serves as the filtering operation. Because A(p) = Sa-{p'^), 
the reduction map A is a decomposable positive map, which 
is generally of formed (p) — Si{p) +S2{p^) with5i,2 being 
two CP maps. 

The last example makes use of an unextendible product ba- 
sis [^, a set of orthonormal product basis {|ai)|/3i)}^i of 
Ti-s (8i Ti-L where S < sL and there is no other pure prod- 
uct state that is orthogonal to this set of basis. If we denote 
P^j:\a,){a,\(^ I A) (A I then p = (1 - - S) repre- 

sents an inseparable states with positive partial transposition. 
If we define 



min(a|(/3|P|a)|/3) 



it can be sure that < e < S/sL [|17|]. Denoting po as a 
normalized density matrix in Hs <E) Hl which has the property 
Tr(po/5) > 0, we define a Hermitian matrix as = P — edpo 
where 

1 = max(a|(/3|po|a)|/3) 

a \a}\l3) 



and 1 < d < sL. Matrix iJ^ is not positive semidefinite since 
TriJe/5 = —edTT{pop) < and for an arbitrary pure product 
state {a\{f3\H^\a)\l3) > 0. If we choose po = h <8) h/sL 
then a non-CP positive map is defined by 



5,(p) = ^T,pi;t-eTrp, 



(13) 



where Ti — \ai) {(3* \ is the half density matrix of the product 
base \ai)\(3i). In Ref. m% pQ is taken as a maximally entangle 
state and d = min(s, iTf Positive map is indecomposable 
because Tg ® Se{p) is not positive semidefinite while Tg (g) 
Sd{p) is positive semidefinite for any decomposable map. 

In conclusion, the concept of the half density matrix was 
studied and its applications to the quantum information are 
discussed in some detail. Based on the half-density-matrix 
representation of a Hermitian linear map, we proved that every 
positive map which is not completely positive is a difference 
of two completely positive maps. A necessary and sufficient 
condition for a non-CP positive map is given, which provides 
a way of constructing such kind of maps. Some examples are 
also presented. Further applications of the half density matrix 
in the quantum information and other fields can be expected 
and the understandings of positive maps provided here may be 
helpful the recognition of the inseparable quantum states and 
to the quantification of the entanglement JT^ ] . 
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